Kinetic model of II-VI(OOl) semiconductor surfaces: Growth rates in atomic layer epitaxy 



CO 

o 

O 

CN 

O 

Q 



o 

CO 



T. Volkmann, M. Ahr, and M. Biehl 

Institut fiir Theoretische Physik und Astrophysik, Sonderforschungsbereich 410, 
Julius-Maxmilians-Universitat Wiirzburg, Am Hubland, D-97074 Wiirzburg, Germany 

(Dated: February 2, 2008) 

We present a zinc-blende lattice gas model of II-VI(OOl) surfaces, which is investigated by means of Kinetic 
Monte Carlo (KMC) simulations. Anisotropic effective interactions between surface metal atoms allow for the 
description of, e.g., the sublimation of CdTe(OOl), including the reconstruction of Cd-terminated surfaces and its 
dependence on the substrate temperature T. Our model also includes Te-dimerization and the potential presence 
of excess Te in a reservoir of weakly bound atoms at the surface. We study the self-regulation of atomic layer 
epitaxy (ALE) and demonstrate how the interplay of the reservoir occupation with the surface kinetics results in 
two different regimes: at high T the growth rate is limited to 0.5 layers per ALE cycle, whereas at low enough 
T each cycle adds a complete layer of CdTe. The transition between the two regimes occurs at a characteristic 
temperature and its dependence on external parameters is studied. Comparing the temperature dependence of 
the ALE growth rate in our model with experimental results for CdTe we find qualitative agreement. 

PACS numbers: 81.10.Aj, 68.55.-a, 68.35.Bs, 81.05.Dz 
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I. INTRODUCTION 

Molecular beam epitaxy (MBE) and its variations continue 
to attract significant interest as a technique for the production 
of, for instance, high quality semiconductor films. In addi- 
tion to its practical relevance, epitaxial growth is highly at- 
tractive from a theoretical point of view. It offers many chal- 
lenging open questions and provides a workshop in which to 
develop novel tools for the modeling and simulation of non- 
equilibrium systems in general. Reviews of the experimental 
techniques and theoretical investigations can be found in, e.g., 
Refs. 1,2. 

II- VI semiconductors, as a promising class of materials, 
have been the subject of many experimental studies, see Ref. 
3 for a recent overview. Among other aspects, competing sur- 
face reconstructions and their interplay with growth or subli- 
mation have been investigated. 4-9 As an example we will con- 
centrate on the CdTe(OOl) surface, in the following, but many 
of the results should apply to other II- VI compounds as well. 
The focus will be on atomic layer epitaxy (ALE), a technique 
which aims at self-regulated layered growth by alternating de- 
position of the elements in ultra high vacuum. Specific prop- 
erties of the material system and the attractive clarity of the 
ALE technique allow for an efficient theoretical modeling of 
the growth scenario. 

In earlier investigations we employed two-dimensional lat- 
tice gas models of the CdTe(OOl) surface 10 ' 11 for the inves- 
tigation of the temperature and flux dependence of surface 
reconstructions. Chemical bonds and the influence of the 
underlying crystal were accounted for by effective pairwise 
anisotropic interactions. Sublimation and ALE growth have 
been studied in the framework of a simple cubic model, 17 
which accounted only for the most essential features of CdTe. 
The correct zinc -blende lattice was first considered in Ref. 16, 
but here we extend the model significantly by taking into ac- 
count the formation of Te-dimers at the surface. More impor- 
tantly, we will consider weakly bound excess Te, which pro- 
vides a reservoir of highly mobile atoms. A similar approach 
with weakly bound arsenic has recently been suggested by 



Itoh et al. 19,20 in order to describe the growth on reconstructed 
GaAs(OOl) and InAs(OOl) surfaces. 

Within the framework of our model we demonstrate how 
the surface reconstruction and limited Cd-coverage restrict 
the ALE growth rate at high temperatures. Our model shows 
how this limitation is overcome if weakly bound excess Te is 
present on the surface at low temperatures. We observe a sud- 
den increase of the growth rate with decreasing temperature, 
in qualitative agreement with experiments. 

Due to the increased number of parameters compared with 
our earlier models a quantitative match with experimental ob- 
servations is hard to achieve at the present stage of the model. 
One possibility to overcome this problem is to combine den- 
sity functional calculations (DFT) with Kinetic Monte Carlo 
(KMC) simulations as was done e.g. by Grosse et al. 21 ' 22 in 
the context of growth on III-V(OOl) semiconductor surfaces. 
However, the available DFT input for II-VI(OOl) surfaces, es- 
pecially CdTe(OOl), does not suffice for a full parametrization 
of our growth model. 

The paper is organized as follows: before introducing the 
model and method in section III, we summarize essential 
properties and experimental findings concerning II-VI(OOl) 
surfaces in the next section. In section IV we discuss the sim- 
ulated growth scenario and present the results. The last sec- 
tion summarizes and gives an outlook on perspective projects 
based on this work. 



II. PROPERTIES OF THE MATERIAL SYSTEM 

A detailed overview of II-VI(OOl) surfaces, with emphasis 
on Te-compounds, can be found in Refs. 4,8,9, for instance. 
Whereas we will mostly refer to CdTe, many of the features 
mentioned here apply to other compounds, e.g. ZnTe or ZnSe. 

CdTe crystallizes in a zinc-blende lattice and consists of 
alternating layers of Cd and Te which are parallel to the (001) 
surface of the crystal. In the bulk, each of these layers forms a 
regular square lattice, cf . Figure 1 . Cd atoms are bound to two 
neighboring Te atoms in the layer below with bonds oriented 
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along the [110] direction, which will be termed y-direction for 
short. Accordingly, Te is always bound to two underlying Cd 
along the [110] or x-direction. 

Under vacuum conditions and in absence of particle depo- 
sition, the CdTe(OOl) surface is terminated by a Cd-layer with 
limited coverage p cd < 1/2. At low temperatures T one finds 
a dominant c(2 x 2) vacancy structure, which corresponds to 
a checkerboard like occupation of the available square lattice 
sites. Above T « 570K the surface is dominated by a local 
(2x1) ordering where the Cd atoms arrange in rows along the 
[110] or ^-direction which alternate with empty rows. Despite 
the fact that significant sublimation sets in even below 570K, 
the temperature driven re-ordering of the Cd-terminated sur- 
face has been discussed in terms of a phase transition in effec- 
tive equilibrium. 5 ' 8 ' 10 ' 11 

Electron counting rules, 12 and density functional theory 
(DFT) calculations 13-15 show that the simultaneous occupa- 
tion of nearest neighbor (NN) sites in the y-direction should 
be extremely unlikely within a terminating Cd layer. Such 
pairs of NN-Cd would lead to an unfavorable local concen- 
tration of positive charges. Note that in both the above men- 
tioned arrangements NN pairs of Cd-atoms in y-direction do 
not occur, indeed. DFT calculations show furthermore that the 
c(2 x 2) ordering has a slightly lower surface energy than the 
competing (2 x 1) structure. 13-15 

The c(2 x 2) ordered Cd-termination can be restored in the 
presence of a stationary Cd flux even at high T. 4 ' 8 On the 
other hand, a steady deposition of Te stabilizes Te-terminated 
surfaces. 4 8 At relatively small fluxes and high temperatures 
a coverage p Tc w 1 is observed, hence there is no vacancy 
structure as for Cd-terminated surfaces. The Te atoms arrange 
in rows of dimers along the [110] -direction at the surface, cor- 
responding to a (2 x 1) reconstruction pattern. 

At high flux and low temperatures Te-coverages p Tc w 3/2 
with a (2 x 1) surface symmetry are observed. It has been 
hypothesized that Te atoms form trimers, 4 8 but the precise 
mechanism of binding the excess Te to the surface remains 
unclear. 9 The evaporation of terminating Te-layers after end- 
ing deposition has been studied in the different temperature 
regimes. 4 ' 8 ' 9 The analysis of the corresponding Arrhenius- 
plots supports the hypothesis of a very weak binding of excess 
Te. 

As we will demonstrate in the following, the influence of 
the excess Te is of particular importance for ALE growth. 
One ALE cycle, i.e. a sequence of one Cd and one Te pulse, 
could add a complete layer of CdTe to the system, in princi- 
ple. However, in a wide range of temperatures, one observes a 
growth rate of at most 1/2 layer per ALE cycle. 6 ' 7 This limita- 
tion reflects the restricted coverage of Cd-terminated surfaces. 
Only for sufficiently low T one finds the expected rate of ap- 
proximately 1 layer per cycle, with a temperature driven, sud- 
den transition from one regime to the other. 6 ' 7 We will explain 
within the framework of our model how the above mentioned 
presence of weakly bound Te atoms at the surface allows to 
overcome the limitation. 
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FIG. 1: Solid-On-Solid representation of the zinc-blende lattice. 
Filled / empty circles represent Cd and Te atoms, respectively, their 
size corresponds to the height above the substrate. The figure shows, 
for instance, a single Cd atom (A) on top of a Te-terminated ter- 
race which extends along the [110] direction. Dotted lines represent 
chemical bonds of Cd to Te atoms in the layer underneath, solid lines 
link Te with underlying Cd. Only bonds between visible atoms are 
displayed. Atoms A and D are examples of mobile particles, whereas 
B and E have more than two chemical bonds and are considered im- 
mobile. Sites C and G are available for Cd-adsorption, whereas F 
can be occupied by a Te atom. Particle A could hop, for example, to 
site C in a diffusion move, or onto G descending the terrace edge. As 
an example for Te diffusion, D could hop to site F. 



III. THE MODEL 
A. Representation of the crystal and energetics 

As a simplification we exclude the incorporation of defects 
or vacancies into the crystal. Hence, the bulk can be rep- 
resented in a Solid-On-Solid (SOS) manner, cf. Figure 1. 
Any atom with more than two chemical binding partners, i.e. 
with at least one NN in the layer above, is considered im- 
mobile. When measuring the surface coverage p cd/Tc with 
Cd/Te atoms we count only those atoms which are mobile. 

Every Cd-Te bond in the system contributes an energy e c < 
0. The crystal structure is further stabilized by an effective 
intra-layerNN interaction et of immobile particles of the same 
species. 

Effective interactions of the mobile particles at the surface 
are chosen to reproduce known properties of the system, e.g. 
surface reconstructions. Consequently, mobile Cd and Te are 
treated on different footings. 

As in the planar models of Refs. 10, 1 1 an infinite repulsion 
excludes NN pairs of mobile Cd atoms in y-direction at the 
surface. The repulsion does not apply to pairs of one mobile 
and one immobile Cd, if the latter has three or four Te part- 
ners. An attractive interaction energy e x < between mobile 
Cd-NN in a>direction and a competing interaction Cd < be- 
tween diagonal neighbors (NNN) favor the (2 x 1) or c(2 x 2) 
arrangement, respectively. As demonstrated in Refs. 10,11, 
the temperature dependence of the Cd reconstruction is repro- 
duced qualitatively, if we set \e x \ < 2 |e<j|. 
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tribute an energy e rn c < 0, whereas the energy e 
assigned to a pair which forms a dimer. Note that dimerization 
introduces additional degrees of freedom to the model without 
changing the topology of the lattice. Pairwise interactions ej d 
of neighboring dimers in x-direction favor their arrangement 
in rows along [110], corresponding to a (2 x 1) surface sym- 
metry. For NN pairs of one mobile and one immobile Te we 
assume the same interaction e b as for pairs of incorporated 
atoms. 

In addition to Te particles at regular SOS lattice sites, we 
include the capture of Te atoms in a weakly bound state, de- 
noted as Te* atoms in the following. The energy contribution 
of each Te* atom is e* where |e*| < 2\e c \. This accounts for 
the fact that the binding of a Te* is weaker than that of an 
isolated mobile atom at a regular lattice site which has two 
Cd-Te bonds. 

Following an earlier hypothesis, 4 8 these weakly bound Te* 
might be associated with atoms which occupy Cd-sites tem- 
porarily and hence form a Te-trimer. 18 But, as already men- 
tioned in section II the precise nature of the Te* is still un- 
known. Here, we abandon the assumption of Te-trimers and 
instead include Te* particles in an effective fashion. 

We introduce a reservoir of Te* with 8* N particles, where 
N is the number of atoms in one complete bulk layer of Te. 
The total reservoir occupation 8* is composed of two parts: 
^* = ®*oc + ®mob- Tne quantity 8* loc accounts for the fraction 
of Te* which are localized on top of mobile Cd and thus neu- 
tralize the infinite repulsion between Cd-NN in y-direction. 18 
This effect will be discussed in more detail in the next sec- 
tion. On the contrary, 8* moh corresponds to the remaining Te* 
atoms which are highly mobile and no assumption is made 
about their precise location and the nature of the binding. In 
general, during our simulations 8* loc <C 8* moh holds for most 
of the time. 

For the calculation of the Te coverage p Tc besides mobile 
Te atoms at regular lattice sites both localized and mobile 
Te* atoms are taken into account as well, i.e. p Tc includes 
8* . According to the experimentally observed range of Te- 
coverages 8 ' 9 the total reservoir occupation thus is restricted to 
values < 8* < 1/2. In our model we will exclude larger 
values by a simple cut-off: if 8* has reached its maximum 
value all processes leading to an increase of 9* are forbidden, 
i.e. their rate is set to zero. 

We furthermore assume that an attractive pairwise interac- 
tion of Te* stabilizes the reservoir. In a mean field fashion we 
represent the latter by an energy contribution N t* h 8* 2 with 
parameter e*. < 0. 

In summary, the total energy of a given configuration can 
be written as 



H = e c n c + e b (n% d +nj e + nl e )+e x n° d + e d n% d 

(1) 



-Te.Te + £ Td n Td + e Td n Td 



+ e m n m 



+ e*8*N + e* b 8* 2 N 



which depends on the number of Cd-Te bonds (n c ), bulk NN 
pairs of Cd (n^ d ) and Te (nj e ), NN pairs of one mobile and 
one immobile Te (nj c ), NN and NNN pairs of mobile Cd at 



the surface (n^ d and n^ d ), and NN pairs of mobile Te not 
forming a dimer (n^ c ). The number of Te-dimers is denoted 
as nj , whereas nj d counts pairs of neighboring Te-dimers 
along the x-direction. The last line in Eq. (1) accounts for the 
contribution of the Te* reservoir which is occupied by 8* N 
atoms. 

Throughout this paper we refer to the following parameter 
set: 



1, e .T 



-1.95, e c = -6.5, e* 



e m — -U.o, e b — — i, e x 



5, e b = -0.8, 
-0.4 (2) 



Note that |e<j| only fixes the energy scale in the following. In 
addition we set ks = \td\ = 1 formally and, hence, measure 
the temperature T on the same scale. Later we will comment 
on a potential quantitative comparison with experiments. 



B. Dynamics and growth 

As we are interested in the non-equilibrium dynamics of the 
system under sublimation or growth conditions, specifying the 
energy of a given configuration is not sufficient. In fact, only 
the energy contribution of a few active layers which contain 
mobile particles will be relevant in the following. The kinet- 
ics of the system is governed by the energy barriers which 
have to be overcome in transitions between the different con- 
figurations. 

Unfortunately, the available experimental data and first 
principles results are not sufficient for a systematic fit of our 
model parameters. In particular, reliable evaluations or esti- 
mates of, e.g., diffusion barriers are not available, apart from 
very few exceptions. Therefore, we do not aim at a precise 
quantitative description of the CdTe(OOl) surface at this stage. 

Nevertheless, essential features can be deduced from ex- 
perimental data and physical insight. Important qualitative 
features of the model turn out to be quite robust against vari- 
ations of the parameters, as long as they comply with some 
essential conditions. In our simulations we will use a parame- 
ter set which is mainly based on our previous investigations of 
related models and their comparison with experimental find- 
ings. 

In our previous studies of simpler models, appropriate 
choices of the parameter set enabled us to reproduce quali- 
tatively the observed temperature dependent structure of Cd- 
terminated surfaces in step-flow sublimation. Furthermore, 
a semi-quantitative match of the re-ordering temperature, 
macroscopic sublimation rates etc. has been achieved, see 
Refs. 16,18. 

The microscopic processes considered in our Kinetic 
Monte Carlo (KMC) dynamics are deposition of Cd and Te 
adatoms, diffusion and desorption of mobile atoms, transi- 
tions of Te atoms between regular lattice sites and the reser- 
voir, and desorption of reservoir atoms. 

Deposition of Cd and Te atoms occurs with fluxes F Cd and 
F Te , respectively. First, a lattice site is chosen randomly. 
Then, a search is performed within an incorporation radius 
Tine = 2 for the site with the lowest height where adsorption 
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TABLE I: Thermally activated processes in the simulation and cor- 
responding values of attempt frequencies z/;, pref actors a; and ac- 
tivation energies Ei used for calculating the Arrhenius rates (Eq. 
(3)) as described in the text. Parameters used in our simulations are 

v = 10 12 s~\ v* = 10 8 s~\ B = 7, B* = 2. The energy change 



AHi for a specific event is g 
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byEq. (1) 




process i 


Vi 


Cti 


Ei 


Cd or Te diffusion 


V a 


1 


max {B , B + AHi} 


Cd diffusion (Te* involved) 


Vo 


u mob 


max {B , B + AHi} 


Cd or Te desorption 


Vo 
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AH, 


Te — > Te* nob 


Vo 


1 - 26* 


max {B , B* + AHi} 


Te - Te* oc 


Vo 


1 - 26* 


max {B , B* + AHi} 


Te^ ob -> Te 


V* 


u mob 


max {B* , B + AHi} 


Te ; * oc -> Te 


V* 


l 


max{B r ,B + AHi} 


Te^ o6 desorption 


V* 


u mob 


AH, 



is possible. 18 If this is, for a Cd atom, a site with a neighbor- 
ing mobile Cd atom in y-direction, adsorption is only possible 
if there is either a localized Te* on top of the Cd neighbor, 
already, or a mobile Te* attaches to the impinging Cd. This 
occurs with probability 9* moh . For a Te atom proper adsorption 
sites include, besides regular Te sites, mobile Cd atoms with 
no localized Te* attached as well as empty Cd sites. In the lat- 
ter two cases, the Te then becomes a mobile reservoir Te* with 
probability 1 - 29*, whereas with probability 29* deposition 
is rejected and the particle evaporates. As a consequence, the 
reservoir occupation is limited to values 9* < 1/2. 

The rates for all thermally activated processes are of the 
Arrhenius form 

ti = vi on exp(-Ei/T), (3) 

where Ei and Vi denote the activation energy and attempt 
frequency for process i, respectively. The prefactor on ac- 
counts for rates which depend explicitly on the reservoir oc- 
cupation (cf. Table I). To calculate the activation energy Ei 
we assume the energy of the transition state H t .i between 
the starting (s) and the final (/) configuration of the sys- 
tem to be equal to the greater of the energies H S i and Hfj 
of both configurations plus an additional barrier: 26 H t ^ — 
ma,x{H s s + B Sti ,Hfs + The values of the barrier 

heights B s ,i and Bfj in general are not the same but de- 
pend on the type of process i. For the activation energies 
Ei = H ti i — H St i one obtains 

E, = max.{B ati ,B fii + AHi} (4) 

for diffusion of mobile atoms and transitions of Te atoms be- 
tween regular lattice sites and reservoir states with the energy 
change AHi = Hfs — H s ^ given by Eq. (1). Diffusion in- 
cludes moves to available NN or NNN sites in the same layer, 
as well as hops across terrace edges, cf. Figure 1 . For all dif- 
fusion steps we set B s s = Bfj = B Q in Eq. (4) with barrier 
B D = 7, similar to Refs. 16,18. On the contrary, the assump- 
tion that Te* atoms in the reservoir are less strongly bound 
implies B s j = B and Bf^ = B* < B D for the transition of 



a Te atom into the reservoir and, conversely, B s ^ — B* and 
Bfs = B D for the opposite process (cf. Table I). 

Desorption of atoms from regular lattice sites or mobile Te* 
requires an activation energy Ei = AHi. We consider only 
desorption of single atoms, neglecting the fact that Te proba- 
bly evaporates from the surface in the form of Te 2 molecules, 
in reality. 

For all diffusion processes and the desorption of mobile 
atoms at regular surface sites we assume an attempt frequency 
Vi = v Q = 10 12 s _1 , whereas a lower frequency Vi = v* is 
used for the transition of Te* atoms to regular Te states and 
the evaporation of mobile Te* . 

This choice is based on the reasonable assumption that Te* 
atoms reside in potential energy minima which are much shal- 
lower than those for Te atoms at regular surface sites. 24 Ac- 
cording to transition state theory the attempt frequency v* 
which corresponds to the harmonic oscillation frequency at 
the local energy minimum will become smaller. 25 A similar 
argument was used by Itoh 20 in order to motivate an attempt 
frequency of weakly bound arsenic as low as 10 8 s _1 . In our 
simulations, we will also set v* = 10 8 s _1 . This choice is fur- 
ther motivated by comparison with experimental observations 
as discussed in the following section. 

The formation or breaking of dimers from NN pairs of mo- 
bile Te atoms is taken into account implicitly. A diffusion 
hop of a Cd atom onto a Te dimer, for instance, requires the 
breaking of the dimer which has to be considered in the energy 
balance AHi of the process. The deposition of a Cd atom will 
always break all underlying Te dimers involved. 

On the other hand, dimer formation is assumed to occur 
instantaneously, whenever it is possible due to the resulting 
energy gain. A similar simplification was used in Ref. 23 for 
the description of Si-dimerization on Si(OOl) surfaces. Hence, 
if any event results in two un-dimerized Te atoms which are 
in the correct configuration for forming a dimer pair, a dimer 
bond will be created immediately. For diffusion and desorp- 
tion events, dimer formation will be accounted for in AH. 

The existence of a weakly bound state for Te atoms is es- 
sential to facilitate growth in our model. Consider a Cd- 
terminated surface with c(2 x 2) reconstruction and maximum 
coverage p cd = 1/2: clearly, deposition of Te at a regular lat- 
tice site would violate the SOS condition. On the other hand, 
the incorporation of additional Cd is impossible due to the in- 
finite repulsion of NN neighbors in y-direction. 

It is plausible to assume, however, that an additional Te 
bound to a Cd would neutralize the repulsion as it provides 
three negatively charged dangling bonds. In our model, mo- 
bile Cd can hop to empty sites neighboring another mobile 
Cd in y-direction, provided there is either already a Te* local- 
ized on top of the neighboring Cd, or a mobile Te* attaches 
to the diffusing particle. This Te* then would become local- 
ized itself. The latter scenario is assumed to happen with a 
probability 9* mob which is reflected in the prefactor on of the 
corresponding rate (cf. Table I). 

Thus, we assume that mobile Te* which are not attached to 
Cd atoms are available instantaneously everywhere at the sur- 
face. Note that after the diffusion step the reservoir occupation 
°* = d *ioc + 9 * mob remains the same but the ratio 9* lo J9* mob has 
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increased. If at any point of the simulation a localized Te* is 
no longer needed because for example the underlying Cd atom 
has evaporated, it will become a mobile Te* again. Through 
the mediation of Te* atoms, Cd atoms at the surface can co- 
alesce into islands in spite of the NN repulsion, provided the 
reservoir of Te* is not empty. 

In our simulations we consider also transitions of Te be- 
tween regular lattice sites and reservoir states as well as the 
desorption of mobile Te* into the vacuum. In general, the 
rates of these processes depend on the reservoir occupation, 
i.e. m 7^ 1 (cf. Table I). For instance a mobile Te at the sur- 
face may become a mobile Te* provided the reservoir is not 
fully occupied. This happens with a rate according to Eq. (3) 
with Ui = (1 — 29*) which is zero if the reservoir density has 
reached its maximum value 1/2. Conversely, the transition of 
a mobile Te* to a regular Te state is proportional to the density 
0, 



mob 



of mobile Te* atoms. Localized Te* atoms which are on 
top of Cd atoms in order to neutralize the infinite repulsion 
may become bound at one of the two neighboring regular Te 
sites, if available, thus reducing the density 9* oc . 

The results presented in the following sections were ob- 
tained in simulations with energy parameters according to Eq. 
(2) and the following barriers and frequencies: 



v = 10 12 s 



= lfjV" 1 , B = 7, B* = 2 (5) 



IV. SIMULATIONS AND RESULTS 
A. The reservoir of excess Te* 

Before presenting results of the ALE simulations we dis- 
cuss the properties of the Te* reservoir in greater detail. In 
our model, the occupation 9* results from the complex inter- 
play of deposition, evaporation and transitions between Te* 
and Te states. However, the basic temperature dependence of 
9* can be understood from considering a flat surface in the 
presence of a steady flux of Te atoms. 

After the Te coverage has reached the value p Tc = 1, only 
deposition directly into the reservoir is possible. If we neglect 
the desorption of regular Te from the surface and transitions of 
the type Te — > Te*, the temporal evolution of 9* is governed 
by the balance of direct deposition and evaporation of Te*. 
The rate of the latter is a function of 9* itself and one obtains 
the following simplifying ODE: 



dO* 
~dt 



— » F le O[l/2-6*]-0*v* exp 



+ 0*1*6*1) 



(6) 

Whereas the Te flux increases 9*, desorption of Te* occurs 
with a frequency proportional to 9* multiplied with the Ar- 
rhenius rate for single particle evaporation from the reservoir. 
Clearly, the filling of the reservoir must saturate at 9* « 1/2. 
However, the precise nature of the saturation process is not 
essential for the following arguments and we choose a simple 
cut-off represented by the Heaviside-function 6 [x] . 

Given e*, e*, and F To , the right hand side in Eq. (6) is pos- 
itive for all 9* G [0, 1/2] if the temperature is low enough. 
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FIG. 2: Analysis of the reservoir dynamics for e* = —5, e*. = —8. 

(a) The temperature T res above which an initially empty reservoir 
of Te* cannot be filled. The curve shows the result as obtained from 
Eq. (6). Additionally, symbols indicate the temperature where the 
growth rate drops from « XML to « Q.5ML CdTe per ALE cycle 
in the KMC simulations (cf. section IV B). 

(b) Sublimation of Te* from an initially filled reservoir for three dif- 
ferent temperatures. Solid lines represent the numerical integration 
of Eq. (6) with F Tc = 0, symbols correspond to simulations of the 
situation in the full model. 



Hence, the reservoir will be filled in this case. Above a charac- 
teristic temperature T res the r.h.s. displays two zeros at fairly 
small values of 9*. As a consequence, the differential equation 
has an attractive fixed point with a stationary value 9* « for 

T > T res . 

Figure 2 (left) shows T res as function of F Tc for the choice 
of parameters e* = -5ande* = -8. There is a rapid increase 
of T res for very low values of F Te which becomes weaker 
and weaker for higher fluxes. With F = 5ML/s we obtain 
T res w 0.39, as an example, and at T = 0.40 the stationary 
occupation is as low as 9* w 0.02. Hence, starting from 9* = 
0, the reservoir cannot be filled significantly for T > T res . 
The behavior predicted from Eq. (6) is in excellent agreement 
with simulations of the situation in the full model. We would 
like to mention that in Ref. 6 a simple evaporation model of 
ALE growth was proposed which produces a flux-dependence 
of the characteristic temperature which is similar to the one in 
Figure 2 (left). 

The simplified Eq. (6) is also useful in investigating the ini- 
tial sublimation of Te* from a surface with maximum cov- 
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FIG. 3: ALE simulations with T = 0.44, snapshots of the surface at (left to right) t = Is, 2s, 3s, 4s corresponding to the end of stages 
I,II,ni,IV in the high temperature growth scenario as described in the text. Dark / light grey spheres correspond to Te / Cd atoms, respectively. 
Darker portions of the surface correspond to lower height above the substrate. Te-dimers are represented by smaller lateral distance of NN 
atoms. 



erage p Tc = 3/2, i.e. 9* = 1/2. Figure 2 (right) displays 
the evolution of 9* for e* = — 5, e*. = -8 and three differ- 
ent temperatures as obtained from the numerical integration 
of (6) with F Te = 0. The results from simulations of the 
corresponding situation in the full kinetic model are again in 
excellent agreement with the simplified description, which ne- 
glects evaporation of regular Te and other processes. 

One obtains a well-defined characteristic time r* for the de- 
crease of 9* from 0.5 to, say, 9* = 0.05. The precise value of 
9* is of little relevance as 9* (t) decreases to zero very rapidly, 
cf. Figure 2 (right panel). The temperature dependence of 
r* is very well described by an Arrhenius law of the form 
r* = t e Ea l T with macroscopic activation energy E a w 8.5 
and prefactor r « 3.8 x 10~ 9 s. In Refs. 4,8 experimental 
data is discussed for the initial sublimation of Te starting from 
coverage 3/2. There, the macroscopic activation energy is re- 
ported to be E a « 0.96 eV, providing a potential reference 
point for attempting a quantitative fit of the model. 

Note that the prefactor r D ~ 10~ 9 s obtained from Fig. 2 
(right) is not too different from the time constant 1/v* = 
10~ 8 s assumed for Te* desorption. This should be expected, 
as the latter is the limiting microscopic process of the sce- 
nario. In a naive attempt to extract the prefactor from an anal- 
ogous Arrhenius plot of experimental data 4 we obtain a value 
of about the same order, i.e. r D ~ 10~ 9 s. On the contrary, the 
value from data for Cd desorption 4 corresponds to frequencies 
which are larger by orders of magnitude. Of course, a reliable 
determination of prefactors from experimental data for rela- 
tively small temperature ranges is questionable if not impos- 
sible (and was not attempted by the authors of Ref. 4). Nev- 
ertheless one might interpret the very rough comparison as 
another qualitative justification for choosing v* much smaller 
than other frequencies in the model. 



B. ALE growth 

In order to gain a qualitative understanding of the ALE 
growth scenario we perform simulations with the parameter 
set specified in Eqs. (2,5). All simulation results presented 
here were obtained from systems with 64 x 64 Cd and Te 



atoms in a bulk layer, averaged over 5 independent runs. 

We model a situation in which Cd and Te are deposited 
in alternate pulses of length 0.9s, each with a constant flux 
F cd = F Te = bML/s. The pulses are separated by a dead 
time of 0.1s, hence, the duration of one complete ALE cycle 
is 2s in total. Note that the characteristic time r* for reser- 
voir sublimation is large compared with the dead time interval 
for all relevant temperatures. In the results discussed here, the 
initial surface was perfectly flat and Te-terminated with cov- 
erage p Tc = 1, 9* = 0. Deposition always started with a Cd 
pulse. 

Illustrating surface snapshots of growing surfaces are 
shown in Figures 3 and 4 for the two different growth regimes 
discussed in the following. Additional illustrations and mpeg- 
movies are available from our web pages. 27 

At fairly high temperatures, growth proceeds according to 
the following scenario in our simulations: 

I) The first Cd pulse adds half a layer of Cd to the system. 
The infinite repulsion of NN in y-direction prevents 
coverages p cd > 1/2. The predominant arrangement 
of Cd atoms is in a local c(2 x 2) ordering, i.e. a checker- 
board pattern. The contribution of the (2x1) row struc- 
ture as a thermal excitation increases with higher T 

II) After the first Te pulse and the dead time, half a layer of 
Te has been incorporated. Diffusion of particles leads 
to the formation of islands upon the Te-terminated sur- 
face. Note that the existence of the Te* -reservoir state 
is crucial for this process, as discussed above. Surface 
Te atoms form dimers, preferentially, and these arrange 
in rows along the a>direction, yielding a predominant 
(2x1) symmetry of the surface. Although incorpora- 
tion of Te occurs to a considerable degree through the 
reservoir state Te*, at the end of the dead time one finds 
9* « 0, as discussed above. 

Ill) The second Cd pulse places mobile atoms on top of ex- 
isting islands and in between. The pre-dominant local 
ordering is the (2 x 1) pattern, i.e. rows of Cd alternating 
with empty rows. This is due to the influence of the is- 
land edges on the reconstruction, see Refs. 16-18 for a 



7 




discussion of an analogous effect in layer by layer sub- 
limation. At the end of the dead time, the islands from 
stage II are still present but now they are Cd-terminated. 

IV) During the second Te pulse, atoms are deposited onto 
the Cd-terminated islands and in between. The particles 
rearrange by means of diffusion, leveling off the islands 
and forming a flat Te-terminated surface. At the end of 
the dead time, the initial state of the system is restored 
with one complete layer of CdTe added to the system. 
In the simulations the filling of gaps is incomplete to a 
degree which depends on the temperature and the length 
of the dead time interval. 

Hence, at temperatures T > T res , our model reproduces ALE 
growth in the same double cycle sequence which was hypoth- 
esized and discussed for CdTe and ZnTe in the literature, see 
Refs. 8,9 and references therein. Figure 5 shows the evolution 
of the film height with time in the ALE process at T — 0.44, 
which is well above the characteristic T res w 0.39 obtained 
for^ To = hMLja. 

The picture changes qualitatively for temperatures T < 
T res . Again, a flat Te-terminated surface with coverage p Tc = 
1, 9* = was prepared, and the deposition started with a Cd 
pulse. The surface structures obtained after the first Cd and 
Te pulse are very much alike the ones of stages I and II of the 
high temperature scenario, since initially the Te* reservoir was 
empty. Note however, that due to a decreased mobility of the 
atoms the islands on the Te terminated surface are less com- 
pact. Since T < T res , at the end of the Te pulse the surface is 
not only Te terminated but also the Te* reservoir is filled, i.e. 
9* « 1/2. 

Thus, from this point on the ALE growth scenario proceeds 
along the following lines in our model: 

A) The Cd pulse places atoms onto the islands and in be- 
tween which is analogous to stage III of the high tem- 
perature scenario. But now Cd particles at the surface 
can coalesce, because Te* from the reservoir are avail- 
able to neutralize the NN repulsion of mobile Cd in y- 
direction. Therefore the gaps between the islands are 
leveled out by half a monolayer of Cd atoms and the 
reservoir which provides also approx. 1/2 monolayer of 
Te. The surface is then covered by another half layer of 



Cd atoms arranging preferentially in a c(2 x 2) vacancy 
structure. Thus, the adsorption of one complete mono- 
layer of Cd and a half monolayer of Te which came 
from the reservoir leads to an increase of the surface 
height by 3/4 of a complete CdTe layer. At the end of 
step A) the reservoir is empty. 

B) The following Te pulse leads to a Te terminated sur- 
face with islands, which is identical to the one described 
above. Again, due to low temperature and short dead 
time, also the Te* reservoir is filled. Hence, another 
half monolayer of Te is adsorbed and in total one com- 
plete layer of CdTe added in only one ALE cycle. At 
the end of step B) the reservoir has reached its maxi- 
mum occupation. 

As an example, Figure 5 shows the evolution of the mean 
surface height (h) with time at T = 0.36, i.e. below T res for 
F Tc = hML/s. Note that for the first Cd and Te pulse, the 
curve is essentially identical with that of the high T regime, 
because in both simulations the reservoir was empty, initially. 
Hence the upper curve of Figure 5 corresponds to the sequence 
of stages A/B only for times t > 2s. The mean height as dis- 
played in Figure 5 increases by 3/4 of a complete CdTe layer 
during the Cd pulse A) because (h) takes into account only 
particles at SOS lattice sites, disregarding the Te* reservoir. 

From simulations for different temperatures we have ob- 
tained the average growth rate per ALE cycle. Figure 6(b) 
shows that, with increasing T, the growth rate drops very 
rapidly in the vicinity of the characteristic T res as obtained 
from the considerations of section IV A. The comparison with 
Figure 6(a) shows qualitative agreement with experimental re- 
sults for an analogous growth scenario. 6 Very similar data has 
been published by Hartmann et aJ. 7 and is also discussed in 
Ref. 8. Note that in both experiments growth rates above 1 
layer of CdTe per ALE cycle were observed at very low tem- 
peratures. This is presumably due to the formation of Cd crys- 
tallites at the surface and cannot be reproduced in our model. 9 

The temperature driven break-down of the growth rate is 
due to the fact that the reservoir cannot be filled significantly 
for T > T res . Simulation results for various values of the flux 
strongly support this assumption as can be seen from Figure 2 
(left): the symbols which indicate the temperature where the 
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FIG. 5: Mean height of the film vs. time in the simulations of ALE 
for two different substrate temperatures. The pulse duration is 0.9s 
with F Tc = F Cd = 5ML/s followed by a dead time of 0.1s. (h) 
takes into account only particles at regular lattice sites, excluding the 
reservoir of Te* . The upper curve evolves according to the low tem- 
perature scenario (T = 0.36) with steps A and B whereas the lower 
corresponds to T = 0.44 with the double cycle sequence I,II,III,IV 
as described in the text. The scenarios essentially coincide for the 
first 2s because in both cases the reservoir is empty before the first 
Cd-pulse is deposited. 



growth rate drops in the simulations lie almost exactly on the 
curve that marks T res . 

In principle, the re-evaporation of Te* during the dead 
time interval could be an alternative cause for the reduced 
growth rate. Note, however, that even for a flux as high as 
F Tc = lOML/s with T res w 0.41, one finds that the charac- 
teristic time for reservoir sublimation is r* « 4s, cf. Figure 
2 (right panel). Hence, it appears safe to say that during dead 
time intervals shorter than, say, Is, 9* will remain close to the 
maximum value 1/2, provided all other conditions allow for 
the filling of the reservoir during a Te pulse. 

As a consequence, we expect that the characteristic value 
of T which marks the drop of the growth rate in Figure 
6 should be essentially independent of the pulse duration 
or dead time interval within a wide range of reasonable values. 

Comparing the experimental results by Hartmann et al. 7 and 
Faschinger and Sitter 6 for the ALE growth rate one finds that 
for the latter the temperature where the growth rate drops to 
approx. 1/2 layers of CdTe per cycle is shifted to higher tem- 
peratures by about 4QK. Within the framework of our model 
this can be explained if one takes into account the different 
experimental conditions. Whereas Hartmann et al. combined 
particle fluxes of about Q.bML/s with 8s pulse time and Is 
dead times, the flux used by Faschinger and Sitter was 2ML / s 
with a pulse duration of 0.8s and a dead time of 0.2s. As 
discussed above, the difference of the dead times can be ne- 
glected, but for the given fluxes one reads off from Fig. 2 
(left) a significant increase of T res by about 9%. As demon- 
strated, T res practically coincides with the temperature where 
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FIG. 6: Growth rates in layers of CdTe per ALE cycle vs. tempera- 
ture T. 

(a) Experimentally observed growth rates of CdTe(OOl) as published 
by Faschinger and Sitter, plot generated after Fig. 3 in Ref. 6. 
The pulse duration was 0.8s with a total deposition of 1.5 layers 
of each element, followed by dead time intervals of length 0.2s. The 
most prominent feature is the sudden decrease of the growth rate 
at T ?s 290°C One additional data point has been omitted here 
(growth rate 1.4 layers per cycle at 230° C). 

(b) Average growth rate in our simulations with parameters as in 
Fig. 5 for different temperatures. T is measured in dimensionless 
units of |e<j| = L The growth rate drops at T w T res w 0.39 for 
F Tc = 5ML/s. 



the growth rate drops which consequently is expected to in- 
crease by the same amount. The difference of the absolute 
temperatures observed in Refs. 7 and 6 is about 8% (of the 
lower one) which is indeed comparable with our results. 

For very high T, sublimation of the crystal dominates over 
the incoming flux and growth becomes impossible. The ALE 
growth rate vanishes at a temperature which is essentially in- 
dependent of the properties of Te* atoms. 



V. SUMMARY AND OUTLOOK 

In summary, we have presented a detailed lattice gas model 
of II- VI semiconductor compounds. In contrast to many sim- 
pler growth models, it takes into account the correct zinc- 
blende lattice as well as essential, material specific effects 
which have been discussed in terms of the CdTe(OOl) surface, 
as an example. The model reproduces, for instance, the re- 
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construction of Cd-terminated surfaces by means of vacancy 
structures and its dependency on temperature and flux. This 
includes the temperature driven transition from the c(2 x 2) re- 
constructed surface to a disordered state with local (2x1) or- 
dering. This re-arrangement of the surface is observed under 
step flow sublimation, as demonstrated earlier in Refs. 16,17, 
and a semi-quantitative agreement with experimental data can 
be achieved. 

Here, in addition, the dimerization of Te at the surface is 
introduced to the model. Thus, the (2x1) reconstruction of 
Te-terminated surfaces is reproduced. We demonstrate fur- 
thermore that the presence of weakly bound Te* atoms is nec- 
essary to facilitate growth by, e.g., atomic layer epitaxy. Al- 
though the precise nature of such a weakly bound Te* state at 
the surface is unknown, its existence is clearly evident from 
various experimental observations. 8 ' 9 

By means of Kinetic Monte Carlo simulations of our model 
we show that the reservoir of Te* also explains the experimen- 
tally observed temperature dependence of the ALE growth 
rates. Only below a characteristic temperature a significant 
amount of Te* is present and facilitates growth rates of ap- 
proximately 1 layer of CdTe per ALE cycle. As in the cor- 
responding experiments, we observe a sudden transition from 
the self-limited growth rate of about 1/2 layer per ALE cycle 
at high T to the low temperature regime with a rate close to 1 
complete layer per cycle. 

Our analysis shows that the characteristic transition tem- 
perature should depend only weakly on the model (or experi- 
mental) parameters within a wide range of reasonable choices. 
The key dependence is on the particle flux during Cd and Te 
pulses, respectively, and our model reproduces its effect on 
the transition qualitatively correct. 



As an attempt towards a more quantitative comparison with 
experiments, one might identify the temperature T res w 0.365 
for a flux of 2ML/s in our model (cf. Fig 2) with the exper- 
imental value T w 290°C by Faschinger et al. 6 Thus, the 
energy scale a may be fixed, and we obtain, for example, a 
value of E a « 1.13 eV for the macroscopic activation en- 
ergy for desorption of Te* atoms which roughly agrees with 
the experimentally determined 0.96 eV. 4,& Furthermore, the 
value of T res for our simulations with 5ML/s translates into 
T res « 329°C which agrees well with results of the simple 
evaporation model from Ref. 6. 

We thus believe, that, in principle, our model offers the po- 
tential for a quantitative match with experimental data. In 
simpler versions of the model it was possible to reproduce 
the temperature dependence of the surface reconstruction as 
well as macroscopic sublimation rates on a semiquantitative 
level. 16 The larger number of parameters in the current exten- 
sion, however, requires further input from experiment. More- 
over, we hope to obtain reliable estimates of microscopic ac- 
tivation barriers from first principle quantum chemical calcu- 
lations as are already available for III- VI semiconductors. 21 ' 22 
Such calculations will be extremely useful and should also 
shed light on the nature of the weakly bound Te* states on the 
surface. 
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